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Abstract 

A star-triangle mapping transformation is used to establish an exact correspondence 
between the spin-1/2 Ising model on the Shastry-Sutherland (orthogonal-dimer) lat- 
tice and respectively, the spin-1/2 Ising model on a bathroom tile (4-8) lattice. Exact 
results for the critical temperature and spontaneous magnetization are obtained and 
compared with corresponding results on the regular Ising lattices. 
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1 Introduction 



Over the last six decades, much effort has been devoted to determine a critical- 
ity and other statistical properties of various lattice-statistical models, which 
would enable a deeper understanding of order-disorder phenomena in solids. 
The planar Ising model takes a prominent place in the equilibrium statistical 
mechanics from this point of view as it represents a rare example of exactly 
solvable lattice-statistical model since Onsager's pioneering work [1] . It is note- 
worthy that early exact studies of the planar Ising model have been performed 
on regular lattices - square [1], triangular and honeycomb [2], having all sites 
as well as all bonds equivalent. Exact calculations performed on these lattices 
revealed the universal critical behavior supporting the universality hypothe- 
sis, in fact, all three regular lattices exhibit a second-order phase transition 
characterized by the same set of critical indices. 
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A remarkable generalization has been achieved when the planar Ising mod- 
els with non-crossing bonds were shown to be, at least in principle, exactly 
solvable within the Pfaffian method [3]. In spite of the universality in their 
critical behavior, it is still of great research interest to calculate accurately 
the critical point and other statistical properties for a specific lattice with 
particular interactions. Exact studies of the Ising model on the Union Jack 
(centered square) lattice proved, for instance, a possible existence of reentrant 
transitions in a presence of competing interactions [4]. It should be mentioned, 
however, that a precise treatment of planar Ising models is usually connected 
with considerable difficulties related to sophisticated mathematical methods 
used. On the other hand, adopting the well-known formalism based on map- 
ping transformation technique [5] largely simplifies this treatment. By the use 
of mapping, exact solutions of the Ising model on a pentagonal lattice [6] and 
two topologically different 4-6 lattices [7]- [8] have recently been derived. 



From a topological viewpoint, the Archimedean (uniform) lattices attract an 
appreciable attention because they represent infinite tessellation in which all 
sites are equivalent [9]. In addition to a subset of three regular lattices, the 
other eight Archimedean lattices being irregular since their tilings consist of 
two or more regular polygons. Although the Ising model on some of these 
'semi-regular' lattices have exactly been treated long ago, e.g. on Kagome 
[10], expanded Kagome (3-12) [11] or bathroom tile (4-8) [12] lattices, some of 
them seem not to have been dealt in the literature so far. The primary goal of 
this work is to find out an accurate solution for the spin-1/2 Ising model 
on the Shastry-Sutherland lattice [13] (to be further abbreviated as SSL) 
schematically shown on right-hand side of Fig. 1. The SSL is five-coordinated 
Archimedean lattice formed by square plaquettes each surrounded by four tri- 
angular plaquettes. Consequently, there are two non-equivalent bonds in the 
SSL which may be a common edge either of two triangles, or a triangle and a 
square. Notice that the lattice topology of SSL can also be viewed as a pla- 
nar network constituted by mutually orthogonal neighboring dimers (see Fig. 
1). It is worthwhile to mention that a theoretical investigation of the Ising 
model on the SSL is interesting from the academic point of view as well as 
in connection with a recent experimental observation of insulating molecular- 
based compounds SrCu 2 (B0 3 ) 2 [14], Nd 2 BaZn0 5 [15] and [Ni(mettn)(N 3 )] [16] 
whose magnetic structure matches a requirement of the SSL topology. 



The organization of this paper is as follows. In the next section we derive 
an exact solution for the spin-1/2 Ising model on the SSL within an exact 
mapping to a corresponding bathroom tile (4-8) lattice. Exact results for the 
critical temperature and spontaneous magnetization are studied particularly 
in Section 3. Finally, some concluding remarks are given in Section 4. 
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2 Model and its exact solution 

Let us begin by considering the anisotropic bathroom tile (4-8) lattice schemat- 
ically illustrated on the left-hand side of Fig. 1. While all sites of the 4-8 lattice 
are equivalent, having coordination number three, the lattice is irregular due 
to the non-equivalency of its bonds. In this respect, two different coupling 
parameters are assigned to the bonds within square plaquettes (R^) and re- 
spectively, to the links between them (Ri). Another important topological 
feature of the 4-8 lattice is that it belongs to 'loose-packed' lattices, i.e. it can 
be divided into two equal sub-lattices A and B in such a way that all nearest 
neighbors of A-sites belong to the sub-lattice B and vice versa. To emphasize 
this topological feature, the A-sites (5-sites) are displayed in Fig. 1 as empty 
(solid) circles. On assumption that all sites are occupied by spin-1/2 atoms, 
the Ising Hamiltonian defined upon the underlying 4-8 lattice reads: 

N 2N 
(ij)cC (ij)cS 

where S" = ±1 (a = A, B) is the Ising spin located at ith. lattice point, its 
superscript specifies the sub-lattice to which it belongs, 2N denotes the total 
number of lattice sites and the summations run over two sets of bonds, the 
bonds within the square plaquettes and the links between them, respectively. 

It is of principal importance that the 4-8 lattice can be divided into two 
interpenetrating sub-lattices. Owing to this fact, a summation over the spins 
of one sub-lattice (say sub-lattice B) can be performed before summing over 
the spins of another sub-lattice (sub-lattice A). Accordingly, the partition 
function of 4-8 lattice can be rewritten in this useful form: 

^4-8 = E fi E exp[5f + L 2 {S£ + Sf)], (2) 

{S A } i=l Sf =-1 

where L\ = Ri/(k-QT) and L 2 = i^A&B^) are effective coupling parameters, 
&b is Boltzmann's constant and T the absolute temperature. Further, the 
product is over all 5-sites and the symbol Z){s^} means a summation over 
all possible spin configurations on the sub-lattice A. Apparently, this form 
justifies performance of the familiar star-triangle mapping transformation [5]: 

E exp[Sf(L 1 S; + L 2 (^ + ^)] 

5f=-l 

= AeM*YS£sf + K 2 Sf(S£ + Sf)}. (3) 
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A physical meaning of the mapping (3) is to remove all interaction parameters 
associated with the B-site located at ith lattice point and to replace them by 
new effective coupling parameters K 1 = Jmtra,/ (k B T) and K 2 = Jmtcr/(k B T) 
between the remaining A-sites. When the mapping (3) is applied to all B- 
sites, the 4-8 lattice is then mapped on the corresponding SSL (Fig. 1). In 
what follows, we will refer to new effective interactions K\ and K 2 of the SSL 
as to the intra- and inter-dimer couplings, respectively. Note that the mapping 
parameters A, K\ and K 2 are unambiguously given by the transformation (3), 
which must be valid for any configuration of Sf, S£, Sf spins. Hence, 



A = 2[cosh(Li + 2L 2 ) cosh(Li - 2L 2 ) cosh 2 (Li)] 1/4 , 



K 1 = -ln 



Ko = -\n 



cosh(Li + 2L 2 ) cosh(Li - 2L 2 ) 



cosh (Li) 



cosh(Li + 2L 2 ) 



cosh(Li — 2L 2 ) 



(4) 
(5) 

(6) 



It is worthy to notice that a substitution of the transformation Eq. (3) into 
the formula (2) establishes a relationship between the partition function of 
the 4-8 lattice and its corresponding partition function of the SSL lattice: 



Z 4 _ 8 (Li,L 2 ) =A"Z SSL (K 1 ,K 2 ). 



(7) 



Of course, similar mapping relations can also be established for other im- 
portant quantities such as Gibbs free energy, internal energy, magnetization, 
correlation functions, specific heat, etc. For instance, with the aid of Eq. (7) 
it can readily be proved that a magnetization of the 4-8 lattice, characterized 
by the effective interactions L\ and L 2 , directly equals to a magnetization of 
the SSL with the coupling parameters K\ and K 2 given by (5)-(6): 



m 4 _ 8 (Li, L 2 ) = m SS L(K 1 , K 2 ). 



(8) 



Now, we will find the inverse transformation in order to express an exact 
solution of the SSL in terms of the known exact solution of the 4-8 lattice [12]. 
After elementary calculations one finds inverse transformation formulae: 



In 



e^i cosh 2K 2 ) 2 - 1 + e^ 1 sinh 2K 2 



e*i cosh 2fT 2 ) 2 - 1 - e K ^ sinh 27^ 



e Kl cosh 2fT 2 ± ^J(e K ^ cosh 2K 2 f 



(9) 



(10) 



It should be stressed, however, that a validity of the mapping relation (9) 
necessitates K 1 > 0, i.e. the mapping between the SSL and 4-8 lattice holds 
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for the ferromagnetic intra-dimer interaction only. On the other hand, a sign 
ambiguity in the mapping relation (10) is of less importance, since basic fea- 
tures of the 4-8 Ising lattice such as critical frontiers, partition function or 
Gibbs free energy remain invariant under the sign change transforming merely 
J 2 — > — L 2 . With regard to this, the partition function and magnetization of 
spin-1/2 Ising model on the SSL lattice can immediately be calculated from 
the corresponding quantities of the 4-8 lattice [12] with the help of Eqs. (7) 
and (8) when exploiting the transformation formulae (9)-(10). 



3 Results and Discussion 

Let us proceed to a discussion of the most interesting numerical results. For 
simplicity, we shall confine our attention to a particular case when both intra- 
as well as inter-dimer coupling constants of the SSL are assumed to be ferro- 
magnetic. It should be mentioned, however, that the mapping relations (9)- 
(10) are responsible merely for a sign change of the coupling parameter L\ 
when transforming K 2 — > —K 2 , what means that the presented results for 
critical boundaries show the critical temperature of SSL with a ferromagnetic 
intra-dimer and an antiferromagnetic inter-dimer interaction as well. 

The critical temperatures of the anisotropic square, SSL and triangular Ising 
lattices are displayed in Fig. 2 as a function of the ratio between relevant in- 
teraction parameters. The interaction parameters J m tra and J m ter are defined 
for all three Ising lattices in Fig. 3 representing them as solid and dashed 
lines, respectively. For better orientation, Fig. 2A (Fig. 2B) illustrates critical 
boundaries of the anisotropic Ising lattices when Jintra is stronger (weaker) 
than Ji n ter- As one can see from Fig. 2A, the critical temperatures gradually 
increase as Jinter/ Jintra strengthens as many as their maximum values cor- 
responding to isotropic lattices are reached. Apart from this rather trivial 
finding it is quite remarkable to ascertain that the critical temperature of five- 
coordinated SSL is lower than the one of four-coordinated square lattice when 
Jntcr/ Jntra < 0.458. This difference is with a high certainty attributable to 
different relative occurrencies of the Ji ntr a- and J in tcr-bonds, which are in the 
ratio 1:1 for the square lattice, whereas they are merely in the ratio 1:4 for 
the SSL. Another interesting fact to observe here is a significant departure be- 
tween the critical temperatures of both regular Ising lattices and respectively, 
the SSL in a weak-J in ter region, where the critical temperature of SSL does 
not seem to follow a steep logarithmic increase observed for regular lattices. 

As far as a strong- Jinter limit is considered (Fig. 2B), an increasing strength 
of the ratio Jmtra/ Jinter is then responsible for almost a linear increase of the 
critical temperature of triangular and SSL lattices. Dependences drawn for 
critical temperatures of these lattices starting both from the critical point of 
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isotropic square lattice and achieve their maximum values when approaching 
the isotropic limit. Nevertheless, the best linear fits of both the curves reveal 
that a linear term determining a slope of critical curve is for the triangular 
lattice roughly twice as large as for the SSL. This could be understood as a 
consequence of two times greater number of Ji ntr a-bonds added to a simple 
square lattice when constructing the anisotropic triangular lattice (see Fig. 
3). Finally, it appears worthwhile to remark that our numerical value for the 
critical point of isotropic SSL (k B T c / J = 2.926261) correctly reproduces a 
critical value reported by Thompson and Wardrop [17] several years ago by 
the use of Vdovichenko's combinatorial approach. 

To complete our analysis, the temperature dependence of spontaneous mag- 
netization is depicted in Fig. 4 for several values of the ratio Jmter/ ^intra- As 
could be expected, the spontaneous magnetization varies very abruptly slightly 
below its critical temperature as ~ Ae 1 ^ 8 , where A is termed the critical am- 
plitude and e = (T c — T)/T c being the fractional deviation of the tempera- 
ture from its critical value. This magnetization-type singularity, which is often 
called also as an algebraic branch point, can serve as a convincing evidence 
that the anisotropic SSL falls into a standard Ising universality class indepen- 
dently of the Jinter/ <4itra strength. 



4 Conclusion 

Exact results for the critical temperature and spontaneous magnetization of 
the spin-1/2 Ising model on the anisotropic SSL are obtained by means of the 
star-triangle transformation. The star-triangle mapping largely simplifies an 
exact treatment of the anisotropic SSL in that it relates these exact results to 
the ones of corresponding 4-8 lattice. Although the investigation of spin-1/2 
Ising model on the SSL has primarily been stimulated by a recent exploration 
of molecular-based compounds having the magnetic structure of SSL topol- 
ogy [14]- [16], its significance should also be viewed in connection with manifold 
applications of the Ising model in seemingly diverse research areas. Actually, 
there are various physical motivations to examine the lattice-statistical Ising 
models which proved their usefulness in the realm of statistical physics as suit- 
able models for investigating a phase separation in liquids, lattice gas models 
capturing liquid-gas transitions, lattice models of alloys, models explaining 
cooperative phenomena in various biological and chemical systems, and so on. 
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Figure captions 

Fig. 1 A schematic representation of the mapping scheme. The spatially anisotropic 
bathroom tile (4-8) lattice, formally divided into two equal sub-lattices A 
(empty circles) and B (solid circles), is mapped onto the anisotropic Shastry- 
Sutherland (orthogonal-dimer) lattice with the help of star-triangle trans- 
formation removing all 5-type sites. 

Fig. 2 Exact critical temperatures of the anisotropic square, SSL and triangular 
Ising lattices in dependence on the ratio Jmter/ •Antra (Fig. 2 A) and Jintra/ Winter 
(Fig. 2B), respectively. 

Fig. 3 A schematic representation of the anisotropic square (Fig. 3A), SSL (Fig. 
3B) and triangular (Fig. 3C) lattices. The coupling Ji ntr a (Winter) is repre- 
sented by solid (dashed) line style. 

Fig. 4 The temperature variation of magnetization normalized per one lattice site 
of the SSL when the ratio Jmter/ Antra changes. Temperature is scaled in the 
dimensionless units. 
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